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Abstract 

New ladder operators are constructed for a rational extension of the harmonic 
oscillator associated with type III Hermite exceptional orthogonal polynomials and 
characterized by an even integer to. The eigenstates of the Hamiltonian separate into 
m + 1 infinite-dimensional unitary irreducible representations of the corresponding 
| polynomial Heisenberg algebra. These ladder operators are used to construct a higher- 

■ order integral of motion for two superintegrable two-dimensional systems separable in 

cartesian coordinates. The polynomial algebras of such systems provide for the first 
time an algebraic derivation of the whole spectrum through their finite-dimensional 
unitary irreducible representations. 
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1 Introduction 



In the vast domain of quantum nonlinear oscillators, those constructed by supersymmetric 
quantum mechanical (SUSYQM) techniques play an important role (see, e.g., Refs. jH El 
[3j HI [5] and references quoted therein). In particular, a model that appeared in the early 
90s [7] and was re-discovered several times later on (see, e.g., Refs. [HI El HD1 E]) is of 
special interest because its eigenstates can be written in terms of exceptional orthogonal 
polynomials (EOP) [121 US]) a field that has received a lot of attention during the last few 

years HH HS1 HH1 H3 HS1 HS1 1201 EH [221 123]- 

Ladder operators for harmonic oscillator supersymmetric partners are usually con- 
structed by combining the oscillator creation and annihilation operators with the super- 
charges [31 HI [8] or combinations of the latter [51 EH EH] . Together with the Hamiltonian, 
such operators close a polynomial Heisenberg algebra (PHA), which may have infinite- 
dimensional, as well as finite-dimensional unitary irreducible representations (unirreps) 
[21 [21] ■ This construction was carried out for the rationally-extended oscillator model 
referred to above and the corresponding PHA was shown to have two unirreps, an infinite- 
dimensional one spanned by all excited states and a one- dimensional unirrep spanned by 
the ground state [8]. 

Apart from their own interest, ladder operators are also very useful in other contexts such 
as nuclear physics, quantum chemistry or condensed matter. They have also applications 
in the context of mathematical physics and more particularly in the field of superintegrable 
higher-dimensional systems. Considering more specifically the case of two-dimensional 
Hamiltonians, when one leaves the realm of well-studied quadratically superintegrable ones, 
i.e., those allowing two second-order integrals of motion (see, e.g., Refs. [26l E7J EEJ EH [30] 
and references quoted therein) , the direct approach for determining the integrals of motion 
becomes more and more difficult as the order of the integrals increases. This is clearly 
shown in recent works on systems with third-order integrals of motion [3H [32], [33l [34] . For 
such a reason, some other approaches, based on ladder operators [35], recurrence relations 
[36] , or SUSYQM [371 EH] , have been proposed. 
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In two recent studies, the recurrence relation [39J and the ladder operator [ID] methods 
have been applied to construct new superintegrable systems connected with EOP families. 
In the latter work, in particular, some two-dimensional systems obtained from the above- 
mentioned rationally-extended oscillator were considered, but the results did not prove 
entirely satisfactory because it was not possible to derive the whole energy spectrum from 
the representations of the polynomial algebra generated by the integrals of motion, in 
contrast with what happened for some other examples. 

The purpose of the present letter is to introduce some novel ladder operators for such 
a rationally-extended harmonic oscillator and to illustrate their usefulness for providing an 
adequate approach to those superintegrable systems. 

2 New ladder operators for a rationally-extended har- 
monic oscillator 

Let 

H^ = -^ + V (± \x)-E m , V^\x) = W 2 (x)TW'(x) + E m , W{x) = -(0 m (x))' (1) 

be a pair of partner Hamiltonians in first-order SUSYQM [41] . where V^ + '(x) = x 2 (— oo < 
x < oo) is the harmonic oscillator potential, while the factorization function and energy 
(such that # (+) 0m = 0) are chosen as 4> m {x) = li m (x) exp(x 2 /2), E m = —2m — 1, with m 
even [TTJ. Here H m ( x) is a pseudo-Hermite polynomial, defined by % m (x) — ( — i) m H m (ix^) 
in terms of a standard Hermite one. The two Hamiltonians intertwine with 

A = ^- + W(x), A^ = -^- + W(x), W{x) = -x-^, (2) 
ax ax H m 

as AH^ +) = H^A, A ] H^ = H^A\ and the partner potential 

V ( -\x) =x 2 -2 

is a rationally-extended harmonic oscillator considered in many works P [TJ El El QUI [HI EE21 
£3]. 
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Since, for even m, T-L m (x) is strictly positive on the whole real line, the partner H^' has 
an extra bound state below the oscillator spectrum, whose eigenfunction is proportional to 
(f>^-(x). As a consequence, the bound-state energies and wavefunctions of the two partners 
are given by [HI 1301 132"] 



Ei +) = 2(z/ + m + l), v = 0,1,2,..., 
£<-) = 2{u + m + 1), v = -m- 1,0,1,2, 



(4) 



and 



if>W(x)=MWH u {x)e-** , v = 0,1,2,... 



1„2 

e 2 X 



= ATH— ^^^(x), !/ = -m - 1, 0, 1, 2, . . . , 

ii-m\X I 



(5) 



with 

J/o fa) = !> y { J+ m +i( x ) = ~'Hm(x)H v+1 (x) -2rnH m -i(x)H u (x), v = 0,1,2,..., (6) 
and 

AC (+) =(v^2^!)- 1/2 , ,/ = 0,1,2,..., 

(om |\ V 2 (7) 
, AA-^I^+V + ^ + lH" 172 , 1/ = 0,1,2,.... 

The polynomials y^ n \x) } n = v + m + 1, v = — m — 1, 0, 1, 2, ... , are known as the (type 
III) Hermite EOP [T2J [T3] . They form an orthogonal and complete set with respect to the 
positive-definite measure exp(— x 2 ) (H m (x)) 2 dx. 

In Ref. [40J, we considered ladder operators = Aa) A\ b = AaA^ for H^' obtained 
from the creation and annihilation operators a* = —d/dx + x, a = d/dx + x, valid for H( + \ 
and the supercharge operators A\ A in the standard way [31 HJ [8]. Here we plan to build 
some alternative ladder operators c', c. 

For such a purpose, let us first show that one can go from to (up to some addi- 
tive constant) by another path using m first-order SUSYQM transformations characterized 
by the supercharges 

A i = -f+W t {x), Al = ~ + Wi(x) } W i (x)=x + '^-^, t = l,2,...,m. (8) 
dx dx Hi-i Hi 
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On denning 
Hi 



d 2 



dx 2 



+ x z 



U n /U 1 



i — 1, 2, . . . , m + 1, 



(9) 



we indeed get A\Ai = Hi and A\aI = H i+1 + 2 for i = 1, 2, . . . , m, which implies that 
= (# i+1 + 2)ii and A\(H i+1 + 2) = i^ij. Since 

if« = ^ + 2m + 4, if ( ^ = # m+ i + 2m + 2 



(10) 



we infer that 



A m ■ ■ ■ A 2 A 1 H {+) = (#H + 2m + 2)A m ■ ■ ■ A 2 A U 

H^A[A\ ■ ■ ■ Al = A\A\ ■ ■ ■ Al(H^ + 2m + 2), 

which proves the above assertion. It is worth observing here that the operators Ai, A\, 
i = 1, 2, . . . , m, and Hi, i = 2, 3, . . . , m, are only auxiliary operators, some of which 
are singular at x — 0. This means that we actually have a dressing chain of Hamiltonians 
[251113]. 

This chain of m first-order SUSYQM transformations from to if' - ) + 2m + 2 can be 
combined with the transformation from H^ to to provide some raising and lowering 
operators for H^', 

c f = AA\A\---Al, c = A m ---A 2 A 1 A\ (12) 

which are (m + l)th-order differential operators. From the set of intertwining relations 
satisfied by A^ and Ai, i = 1, 2, . . . , m, it is indeed easy to prove that cH^ = (if + 
2m + 2)c or 

(13) 




The operators H^ \ c\ and c fulfil the commutation relations 

[H { -\ c f ] = (2m + 2)c\ [H<-\ c] = -(2m + 2)c, 

[c, c f ] = Q(H^ + 2m + 2) - Q(# M ), 

where 

m 

Q(# (_) ) = H { - ] - 2m - 2 - 2z) 

is a (m + l)th-order polynomial in H^~\ They therefore form a PHA of mth order [3, 



(14) 



(15) 
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Fig. 1. Energy spectrum of H^~' and action of on the eigenstates for m = 2. The v 

values are indicated on the right. 

The action of the raising operator on the eigenfunctions ipl~\x) of H^~\ defined in 
©-(Cj), can be easily calculated and is given by 

ctyii-i= [2™ +1 (m+l)!] 1/2 4 _) , 

ctyH = _ [2 m +> + m)( v + m - 1) • • • (u + l){u + 2m + 2)] 1/2 ^i+i, (16) 
i/ = 0,1,2,.... 
For that of c, we get 

= 0, i/= -m- 1,1,2, ...,m, (17) 

as well as the Hermitian conjugate of Eq. (TIB]) . We conclude that the PHA generated by 
H^~\ c\ and c has m + 1 infinite-dimensional unirreps spanned by the states {' ? /'j-+( m+1 )j I 
j = 0, 1, 2, . . .} with i = —m — 1, 1, 2, . . . , m, respectively. 

The spectrum of H^~' and the action of on the eigenstates is displayed in Fig. 1 for 
the m = 2 example, in which case the PHA has three infinite-dimensional unirreps. 
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3 Application of the new ladder operators to some 
superintegrable two-dimensional systems 



Let us consider some two-dimensional Hamiltonians allowing separation of variables in 
cartesian coordinates, 

H = H x + H y = - ^ + V x (x) + V y (y), (18) 

and let us assume that there exist ladder operators (a x , a x ) and (a y , a y ) in both axes that 
are differential operators of order k\ and k 2 , respectively, and satisfy the defining relations 
of two PHA's, 

[H x , a x ) = X x a\, [H x , a x ] = -X x a x , [a x , a x ] = Q(H X + X x ) - Q{H X ), 

(19) 

[H y , a y ] = AyaJ, [Hy, a y ) = -X y a y , [a y , a y ) = S(H y + X y ) - S(H y ). 
Here X x and X y are constants, while Q(H X ) and S(H y ) are polynomials. 

The separation of variables in cartesian coordinates implies the existence of a second- 
order integral of motion H x — H y , showing that the two-dimensional system (f!8l) is in- 
tegrable. From the ladder operators, one can construct additional polynomial opera- 
tors commuting with H, and a r ^}a} y ' 1 , where ri\, n 2 G Z + are chosen such that 
n\X x = n 2 X y = X [35]. Hence system (fl8l) possesses three algebraically independent inte- 
grals of motion and is superintegrable. It is worth stressing that from ladder operators of 
rather low order k\, k 2 , one generates integrals of motion of higher order k\in\ + k 2 n 2 in a 
nice factorized form that would be difficult to obtain in a direct approach. 

The integrals of motion 

K=^{H X -H y ), I + = a^a n y -, J_ = o^af" (20) 
generate the polynomial algebra of the system 



[K, J±] = ±J ± , [/_,/+] = F num (K + 1,H)- F nijm (K, H), 

F ni:m (K, H) = fj Q (^ + XK - (n x - i)X x ) fl S (^ - XK + jX y \ 

i=l ^ ' 9=1 ^ ' 



(21) 



r- 

which is of order k\n\ + k 2 n 2 — 1. Such a polynomial algebra is a deformed u(2) algebra 
and its finite-dimensional unirreps can be found by realizing it as a generalized deformed 



oscillator algebra {6*, 6, N}. The operators 6* = I + , b — I-, N — K — u and Q(H, u, N) 
F nit n 2 (K, H) indeed satisfy the defining relations of such an algebra 



[N,b t }=b\ [N,b] = -b, 6*6 = <$>(H,u,N), bb l = $(H, u, N + 1), (22) 

where u is some constant and &(H, u, N) is called "structure function" . If the latter satisfies 
the properties 

$(E,u,Q)=Q, $(E,u,p+ 1) = 0, $>{E,u,n)>0, n = l,2,...,p, (23) 

then the deformed oscillator algebra has an energy- dependent Fock space of dimension p + 1 
with a Fock basis \E, n), n — 0, 1, . . . , p, fulfilling 

H\E,n) = E\E,n), N\E,n) — n\E,ri), b\E,0) = 0, b*\E,p)=0, 

(24) 

tf\E,n) = y / <^(E,u,n + l)\E,n + 1), b\E,n) = y/$(E, u, n)\E, n - 1). 
These relations can be used to obtain the (p + l)-dimensional unirreps of the polynomial 
algebra (|2"T1) and the corresponding degenerate energy spectrum of the system. 

Let illustrate this method on two systems already considered in Ref. [ID] by taking 
advantage this time of the new ladder operators introduced in Sec. 2. 

3.1 Combination of a rationally- extended oscillator with a stan- 
dard one 

Let us consider the two-dimensional system given by ([18]) with respectively in the x-axis the 
superpartner of the harmonic oscillator considered in Sec. 2 and in the y-axis the harmonic 
oscillator itself, 



d 2 

H x = H(-> -2m - 1 = — — -+x 2 -2 

ax 2 

d2 2 



-X/" /'J-/' x 2 



Urn V 



1 



(25) 



On taking for (a\.,a x ) the operators (c^,c), defined in (fl2l) . and for (a y ,a y ) standard har- 
monic oscillator creation and annihilation operators, it is obvious that Eq. ( fl9l) is satisfied 
with X x = 2m + 2, X y = 2, and 

m 

Q(H X ) = (H x + 2m + 1) H(H X - 1 - 2i), S(H y ) = H y - 1. (26) 

8=1 



The structure function $(E,u,x) is obtained form Eqs. (121]) (with n\ — 1, n<i — m + 1, 
A = 2(m+ 1)) and ([26]) as 

®{E,u, x) 

= (f + 2 (m + !)0 + u ) + 2 ™ + II (f + 2 ( m + !)(« + «)-!- 2^ 



(27) 



m+l s 

II l K j-2(m+l)(x + u) + 2j-lj . 



From this structure function and the first constraint of Eq. (123]) . we get three different 
types of solutions for the parameter u, 

E 

2(m + = - — - 2m - 1, 

2(m+l)« 2 = -- + 21 + 1, I G {1,2,..., m}, (28) 

2{m + l)u 3 = j + 21-1, I e {1, 2, . . . ,m + 1}. 

The finite-dimensional unirreps are calculated from the two other constraints of Eq. ( 123]) 
and come from the first two types of solutions u\ and U2- They are associated with the 
energies 

E x = 2[(m + l)p+ 1 - k] (29) 

and 

E 2 = 2[(m + l)(p + l) + l-k + l], (30) 

with p G N, k G {1, 2, . . . , m + 1}, / G {1,2, ...,m}, and with appropriate structure 
functions that we skip for brevity's sake. We conclude that there are altogether (m + l) 2 
unirreps characterized by the same peN. 

Let us now show that, in contrast with the previous approach [ID], the present one 
provides all the levels of the physical energy spectrum with their corresponding degeneracy. 
From Eqs. @) and (125]) . the energy spectrum of H is indeed obtained as 

E = E x + E y = 2(y x + u y + l), v x = -m - 1, 0, 1, 2, . . . , ^ = 0,1,2,.... (31) 

Taking iV = v x + v y + 1 and En = 2N, we obtain deg^E^) — 1 il N — —m, —m + 1, . . . , — 1 
and deg(-EAr) = N+l if AT = 0, 1, 2, . . .. On setting v x = (m+l)n x +ax, v y = (m+l)n y +a2, 



with n x , n y G N, a\ G { — m — 1,1,2,..., m}, and 02 G {0, 1, ... , m}, -E can be rewritten as 
E = 2[(m+l)(n x +n y )+ai + a2 + l}. It is then straightforward to see that E\ and E2, defined 
in (12^1) and (13U1) . correspond to with fij. + n 9 = j) 6 N, 02 = m + 1 - G {0, 1, . . . , m}, 
and a>i = — m — 1 or a\ — I G {1, 2, . . . , m}, respectively. 

In the simplest m = 2 case, i.e., for Potential 1 in Ref. |33] , which is one of Gravel's 
systems [32], the polynomial algebra ([2Tj) has nine unirreps for each p G N, associated with 
the energies 6p — 4, 6p — 2, 6p, 6p + 4, (6p + 6) 2 , (6p + 8) 2 , and 6p+10, respectively. On the 
other hand, the sequence of energy levels with their degeneracy is —4, —2, 0, 2 2 , 4 3 , 6 4 , 8 5 , 
10 6 , .... Only the lowest ones belong to a single unirrep (for instance, —4 is obtained from 
6p — 4 with p = and p + 1 = 1), whereas the remaining ones belong to several unirreps 
(for instance, 10 6 is obtained from 6p — 2 with p = 2 and p+1 = 3, 6p + 4 with p = 1 and 
p + 1 = 2, and 6p + 10 with p = and p + 1 = 1). 

In the general case, i.e., for an arbitrary even value of m, a detailed analysis from the 
two solutions Ei and E 2 enabled us to recover the degeneracies. Using N = (m + 1)A + /i, 
we found one unirrep for A = —1 and /i = 1, 2, . . . , m or A = /i = with deg(-Ejv) = 1, 
fi unirreps for A = and fi = 1,2, ... , m with deg(i?7v) = yU + 1 and m + 1 unirreps for 
A = 1, 2, . . . and \i = 0, 1, . . . , m with deg^E^) = (m + 1)A + /i + 1. 

3.2 Combination of two rationally-extended oscillators 

Let us now consider the case where 



with mi and 1112 even and such that m\ > ni2. This includes another system obtained by 
Gravel [32] (Potential 6 in Ref. [33]) for mi = m 2 = 2. 

On taking for (aj., a x ) and (a^, a y ) operators of type (TT2|) (with m — i- m! and m — )■ m 2 , 
respectively) and on noting that Eq. (fl9|) is satisfied with X x = 2mi + 2, A^ = 2m 2 + 2, 
ni = m 2 + 1, n 2 = mi + 1, A = 2(mi + l)(m 2 + 1), and with both Q(H X ) and S(H y ) 




(32) 
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assuming a form similar to Q(H X ) in (1261) (with m->mi and m — » m 2 , respectively), we 
arrive at the following structure function: 



— + 2{m t + l)(m 2 + l)(x + u)- {2m x + 2)(m 2 + 1 - i) + 2m! + 1 



TO2 + 1 

n 

i=l L 

x II ( ^ + 2 ( mi + + l )( x + u ) - ( 2m i + 2)(m 2 + 1 - i) - 1 - 2k 



x 



fe=i 

mi+l 

n 



2~ 



(33) 



— - 2(m x + l)(m 2 + l)(x + u) + (2m 2 + 2)j + 2m 2 + 1 



J] -- 2(mi + l)(m 2 + l)(x + u) + (2m 2 + 2)j - 1 - 2/ 



On proceeding as in the previous case, it can be easily shown that among the four 
different types of solutions for the parameter u, only the first two 



2( mi + l)(m 2 + l)ui = -— + (2mi + 2)(m 2 + 1 - q) - 2m x - 1, 

E 

2{m x + l)(m 2 + l)u 2 = -— + (2m x + 2)(m 2 + 1 - g) + 2r + 1, 



(34) 



(35) 



with g G {1,2,..., m 2 + 1} and r G {1, 2, . . . , mi}, lead to appropriate finite-dimensional 
unirreps of the polynomial algebra ( 12~TT) . The resulting energies are 

E n = 2[(mi + l)(m 2 + l)(p + 2) - {m x + l)g - (m 2 + l)s - (mi + m 2 + 1)], 

#12 = 2 [(mi + l)(m 2 + l)(p + 2) - (mi + l)g - (m 2 + l)s - m x + 1], 
E 21 = 2[(mi + l)(m 2 + l)(p + 2) - (m x + l)g - (m 2 + l)s - m 2 + r], 
E 22 = 2[(mi + l)(m 2 + l)(p + 2) - (mi + l)g - (m 2 + l)s + r + t + 1], 
with p G N, g G {1,2, . . . ,m 2 + 1}, r G {1,2, . . .,m x }, s G {1,2, . . . ,m x + 1}, and t G 
{1,2, .. . ,m 2 }. This time, there are altogether (m x + l) 2 (m 2 + l) 2 unirreps characterized 
by the same p G N. 

On the other hand, from Sec. 2 we know that the physical energy spectrum is given by 



E = E X + E y = 2{y x + v y + 1), v x = -m x -1,0,1,2, 



v = _ m2 _i ; o, 1,2, 



(36) 



On setting v x = (m x + l)[(m 2 + l)n x + a 3 ] + a x , v y = (m 2 + l)[(mi + l)n y + a 4 ] + a 2 , with 
n x , n y G N, a x G {-mi — 1,1,2,..., mi}, a 2 G {-m 2 - 1, 1, 2, ... , m 2 }, a 3 G {0, 1, . . . , m 2 }, 
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04 G {0, 1, ... , mi}, we can rewrite E as E = 2[(mi+l)(m 2 +l)(n a +%) + (mi+l)a3+(m 2 + 
1)04+01+02+1]. We then immediately see that En, E\ 2) E 2 i, and E 22 correspond to E with 
n x +n y = p G N, a 3 = m 2 + l — q G {0, 1, ... , m 2 }, a 4 = mi + 1 — s G {0, 1, . . . , mi}, and ai = 
— mi — 1, a 2 = — m 2 — 1, or % = — mi — 1, a 2 = t G {1, 2, . . . , m 2 }, or a% = r G {1, 2, . . . , mi}, 
a 2 = — m 2 — 1, or ai = r G {1, 2, . . . , mi}, a 2 = i G {1, 2, . . . , m 2 }, respectively. The 
polynomial algebra of the system therefore provides the whole energy spectrum. We have 
checked on several examples that it also accounts for the level degeneracies, which are in 
general obtained through the use of several unirreps. 

Here we have taken the convenient and uniform choice n\ = m 2 + 1 and n 2 = mi + 1. 
However it is worth observing that whenever m x + 1 and m 2 + 1 have a common factor, 
i.e., mi + 1 = \W\ and m 2 + 1 = [iv 2l there exists a simpler choice for n\ and n 2 , namely 
rii = v 2 and n 2 — V\, which would lead to a lower-order polynomial algebra. In any case, it 
is well known that if ladder operators provide an easy method for constructing integrals of 
motion, the resulting algebraic structures are not necessarily the simplest ones that can be 
obtained [5B] . 

4 Conclusion 

In the present work, the construction of ladder operators for a well-known rational extension 
of the harmonic oscillator, associated with type III Hermite EOP and characterized by an 
even integer m, has been reconsidered. Novel operators closing a PHA of mth order have 
been built and it has been shown that the eigenstates of this rational extension separate 
into m + 1 infinite-dimensional unirreps of the PHA. 

Such ladder operators have then been applied to construct a higher-order integral of 
motion for two superintegrable two-dimensional systems separable in cartesian coordinates. 
It has been proved that the polynomial algebras of these systems provide an algebraic 
derivation of the whole energy spectrum through their (p + l)-dimensional unirreps. The 
degeneracy of the energy levels in general results from the union of several unirreps. 

In conclusion, we have shown that as it was already the case for superintegrable systems 
connected with type I or II EOP, a full algebraic treatment may also be found for those 
related to type III ones provided some appropriate ladder operators are constructed. 
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The integrals we constructed with these new ladder operators are of higher order than 
the ones generated by standard ladder operators. These results point out that beyond 
quadratically superintegrable systems the lowest-order integrals do not necessarily provide 
the whole energy spectrum with its degeneracies and that to get the latter one might need 
integrals of higher order. 

In a future work, we hope to be able to carry out a similar study for the systems [10] 
built from extended radial oscillators [18, 19], as well as for those that might be built from 
other extended potentials [UJ H6] . 
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